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The interactions between holes in the Hubbard model, in the low density, intermediate to strong 
coupling limit, are investigated. Dressed spin polarons in neighboring sites have an increased kinetic 
energy and an enhanced hopping rate. Both effects are of the order of the hopping integral and 
lead to an effective attraction at intermediate couplings. Our results are derived by systematically 
improving mean field calculations. The method can also be used to derive known properties of 
isolated spin polarons. 

PACS number(s): 74.20.-z, 02.70.Lq, 71.10.Fd 

The nature of the low energy excitations in the Hub- 
bard model has attracted a great deal of attention. Close 
to half filling, a large amount of work suggests the exis- 
tence of spin polarons, made of dressed holes, which prop- 
agate within a given sublattice with kinetic energy of the 
order of J = ^y- , where t is the hopping integral and 
U the on site Coulomb repulsion. These results are con- 
sistent with similar calculations in the strong coupling, 
low doping limit of the Hubbard model, the t — J model 
|H-H|. There is also evidence for an effective attraction 
between these spin polarons 

In the present work, we analyze the dynamics of spin 
polarons and the interactions between them by means of 
a systematic expansion around mean field calculations of 
the Hubbard model on the bipartite square lattice. Two 
spin polarons in neighboring sites experience an increase 
in their internal kinetic energy, due to the overlap of the 
charge cloud. This repulsion is of the order of t. In ad- 
dition, a polaron reduces the obstacles for the diffussion 
of another, leading to an assisted hopping term which is 
also of the order of t. The combination of these effects is 
an attractive interaction at intermediate values of U jt. 

Use of the Unrestricted Hartree Fock (UHF) approx- 
imation in finite clusters provides a first order approx- 
imation to the spin polaron near half filling. As dis- 
cussed elsewhere, this approximation describes well the 
undoped, insulating state at half filling [fLl|| . A realistic 
picture of the spin wave excitations is obtained by adding 
harmonic fluctuations by means of the time dependent 
Hartree Fock approximation (RPA) [jL2|. At intermedi- 
ate and large values of U/t, the most stable HF solution 
with a single hole is a spin polaron fllj . This solution is 
replaced by a fully ferromagnetic one at sufficiently large 
values of U/t jn|. Approximately half of the charge of 
the hole is located at a given site in the spin polaron so- 
lution. The spin at that site is small and it is reversed 
with respect to the antiferromagnetic background. The 
remaining charge is concentrated in the four neighbor- 
ing sites. A number of alternative derivations lead to a 



similar picture of this small spin bag |fl4|-[l7[ . A similar 
solution is expected to exist in the t — J model. 

As usual in mean field theories, the spin polaron so- 
lution described above breaks symmetries which must 
be restored by quantum fluctuations. In particular, it 
breaks spin symmetry and translational invariance. Spin 
isotropy must exist in finite clusters. However, it is spon- 
taneously broken in the thermodynamic limit, due to the 
presence of the antiferromagnetic background. Hence, we 
do not expect that the lack of spin invariance is a serious 
drawback of the Hartree Fock solutions. In any case, spin 
isotropy can be restored, starting from the mean field 
wavefunction, by projecting out the components which 
do not have a predetermined spin. Results obtained for 
small clusters [ p^!9[ | show a slight improvement of the 
energy, which goes to zero as the cluster size is increased. 
On the other hand, translational invariance is expected 
to be present in the exact solution of clusters of any size. 
Thus, we have improved the mean field results by hy- 
bridizing a given spin polaron solution with all wave- 
functions obtained from it by lattice translations. This 
procedure is equivalent to the configuration interaction 
(CI) method used in quantum chemistry. If the initial 
mean field solution is considered as the "classical" zeroth 
order approximation to the exact solution, this scheme 
can be described as the inclusion of instanton effects, in 
which the spin polarons tunnel between equivalent config- 
urations. Finally, in addition to the previous corrections, 
we can add the zero point fluctuations around the RPA 
ground state |]12fl . This calculation does not change ap- 
preciably the results, although it is necessary to describe 
correctly the long ranged magnon cloud around the spin 
polaron pfJ. 

A schematic picture of the initial one hole and two 
holes Hartree Fock wavefunctions used in this work is 
shown in Fig. [l]. They represent the solutions observed 
at large values of U/t for the isolated polaron and two 
spin polarons on neighboring sites. Actually, charge lo- 
calization is associated to the existence of bound states 
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which split from the top of the lower Hubbard band. 




FIG. 1. a) Sketch of one of the bipolaron solutions, at 
large values of U/t, considered in the text. Circles denote the 
local charge, measured from half filling, and arrows denote the 
spins. There are two localized states marked by the dashed 
line. For comparison, the single polaron solution is shown in 
b). 

Spin polaron wavefunctions localized at different sites 
are not orthogonal. Both wavefunctions overlap and non- 
diagonal matrix elements of the Hamiltonian need to be 
taken into account when mixing between configurations 
is considered. The procedure has been described in de- 
tail in along with the energy improvements (with 
respect to UHF) introduced by this CI scheme. Calcula- 
tions have been carried out on L x L clusters with periodic 
boundary conditions (L < 12) and U > 8t Although 
larger clusters can be easily reached, no improvement of 
the results is achieved due to the short-range charac- 
ter of the interactions. The calculated dispersion of a 
single polaron is shown in Fig. || Because of the antifer- 
romagnetic background, the band has twice the lattice 
periodicity. Exact calculations in finite clusters do not 
show this periodicity, as the solutions have a well defined 
spin and mix different background textures. As clus- 
ter sizes are increased, however, exact solutions tend to 
show the extra periodicity of our results. We interpret it 
as a manifestation that spin invariance is broken in the 
thermodynamic limit, because of the antiferromagnetic 
background. Hence, the lack of this symmetry in our 
calculations should not induce spurious effects. The only 
overlaps and matrix elements which are not negligible are 
those between polaron wavefunctions located in the same 
sublattice. Fig. ^ shows the polaron bandwidth as a func- 
tion of U. It behaves as t 2 /U (the fitted law is given in the 
caption of Fig. ^|). Our scheme allows a straightforward 
explanation of this scaling. Without reversing the spin of 
the whole background, the polaron can only hop within 
a given sublattice. This implies an intermediate virtual 
hop into a site with an almost fully localized electron of 
the opposite spin. The amplitude of finding a reversed 
spin in this new site decays as t 2 /U at large U. We find 
that the polaron band can be very well fitted by the ex- 



pression: efc — eo + ttu cos(k x ) cos(fc y ) + 2£ 2 o[cos(2fc. r ) + 
cos(2k y )} +At22 cos(2fca;) cos(2fc y )-r-4f 3 i[cos(3fc a; ) cos(fc y ) + 
cos(A: x )cos(3fcj,)]. For U = 8t, we get t xx = 0.1899i , 
t 20 = 0.08 73i, t 22 = -0.0 1 36t, and t 31 = -0.0087i. All 
hopping integrals vanish as t 2 /U in the large U limit for 
the reason given above. Also the energy gain with re- 
spect to UHF behaves in this way. Let us mention, 
that all the features reported here are in good agreement 
with known results p]-pj] for both the Hubbard and the 
t — J models. 
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FIG. 2. Quasiparticle band structure for a single hole on 
12 x 12 clusters of the square lattice with periodic boundary 
conditions and U = 8t (filled circles). The continuous line 
corresponds to the fitted dispersion relation (see text). The 
inset shows the bandwidth as a function of t 2 /U; the fitted 
straight line is -0.0224 + ll.llt 2 /f7. 

We now consider solutions with two spin polarons. The 
relevant UHF solutions are those with S z = 0. In order 
to get some coupling, the centers of the two spin polarons 
must be located in different sublattices. The mean field 
energy increases as the two polarons are brought closer, 
although, for intermediate and large values of U, a lo- 
cally stable Hartree Fock solution can be found with two 
polarons at arbitrary distances. We have not attempted 
to do a full CI analysis of all possible combinations of 
two holes in a finite cluster. Instead, we have chosen a 
given mean field solution and hybridized it with all oth- 
ers obtained by a lattice translation or rotation. Clusters 
of sizes up to 10 x 10 were studied which, as in the case 
of the polaron, are large enough due to the short-range 
interactions between different configurations. The ba- 
sis used for the two polarons at the shortest distance is 
shown in Fig. ||. This procedure leads to a set of bands 
for the two hole configurations. The number of bands is 
two or four, depending on the number of different config- 
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urations which can be obtained by rotations at a given 
site. Like in the single polaron case, we obtain a gain in 
energy, due to the derealization of the pair. 
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FIG. 3. Sketch of the bipolaron UHF wavefunctions used 
in this work. Note that the four wavefunctions are obtained 
by successive rotations of tt/2. The complete basis set is pro- 
duced by translation of these wavefunctions through the whole 
cluster. 

The hole-hole interaction, i.e., the difference between 
the energy of a state built up by all configurations with 
the two holes at the shortest distance (separated by a vec- 
tor of the set {1,0}) and the energy of the state having 
the holes at the largest distance possible at a given clus- 
ter is depicted in Fig. ^. Two holes bind for intermediate 
values of U . This happens because the derealization 
energy tends to be higher than the repulsive contribution 
obtained within mean field. The local character of the in- 
teractions is illustrated by the almost null dependence of 
the results shown in Fig. |I] on the cluster size. The en- 
ergy gain of the two holes (with respect to UHF) in the 
two limiting configurations (at the shortest or largest dis- 
tance possible) is given in the inset of Fig. ^. Note that, 
whereas in the case of the holes at the largest distance, 
the gain goes to zero in the large U limit, as for the iso- 
lated polaron, when the holes are separated by a {1,0} 
vector the gain goes to a finite value. This result is not 
surprising, as the arguments given below suggest, and is 
in line with the results for the width of the quasiparti- 
cle band. The numerical results for L—6, 8 and 10 and 
U in the range 8t ~ 5000i can be fitted by the following 
straight lines, 3.965£ + 14A7t 2 /U (holes at the shortest 
distance) and -0.007i + I0.lt 2 /U (holes at the largest 
distance). Thus, total bandwidth of the two bands ob- 
tained for holes in neighboring sites does not vanish in 
the infinite U limit (as the energy gain reported in Fig. 
4) . The internal consistency of our calculations is shown 
comparing the large U behavior of the two holes at the 
largest distance possible with the corresponding results 
obtained for the isolated polaron (compare this fitting 
with that given in the captions of Fig. 0) . 
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FIG. 4. Comparison of the hole-hole interaction (see main 
text for the definition) obtained within UHF (empty symbols) 
and CI (filled symbols) approximations. Results correspond 
to 6 x 6, (circles) 8x8 (squares) and 10 x 10 (diamonds) 
clusters with periodic boundary conditions. The inset shows 
the energy gain due to the inclusion of correlation effects via 
CI for both the configuration of holes located in neighbor- 
ing positions (filled simbols) and holes that are maximally 
separated in the finite size cluster (empty simbols). The fit- 
ted straight lines are: 0.495t + 1.53t 2 /U (filled symbols) and 
-0.002* + 3.78t 2 /U (empty symbols). 

Both the behavior of the quasiparticlc bandwidth and 
that of the energy gain with respect to UHF of the two 
holes at the shortest distance can be understood from 
arguments similar to those used for the single polaron. 
The hopping terms that are proportional to t at large 
U describe the rotation of a pair around the position of 
one of the two holes. Each hole is spread between four 
sites. In order for a rotation to take place, one hole has 
to jump from one of these sites into one of the rotated 
positions. This process can always take place without a 
hole moving into a almost fully polarized site with the 
wrong spin. In the single polaron discussed be- 

fore, the motion of a hole involves the inversion of, at 
least, one almost fully spin polarized site, in the large 
U limit. As a consequence, the derealization of polaron 
pairs on neighboring sites leads to a finite gain in energy, 
even as U — > oo, as opposed to the single polaron case 
p8| or polaron pairs at the largest distance (see inset of 
Fig. @). 

The possibility of hole assisted hopping was discussed 
in j24|, in a different context. It always leads to super- 
conductivity. In our case, we find a contribution, in the 
large U limit, of the type: 

+ i-\,j;sCi+i,y,s + h.c. + perm) (1) 
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This 



term 



admits 



the 



BCS 



decou- 
It favors 



pling At(c[ J . s cl +l j . s )c lJ . s c l . 3+hs + h.c. + . 
superconductivity with either s or d wave symmetry, de- 
pending on the sign of At. Since we find At > 0, d wave 
symmetry follows. 

An interesting question is whether the holes would tend 
to segregate when more holes are added to the cluster. 
In order to investigate this question, we have calculated 
the total energies for four holes centered on a square and 
two separated bipolarons with holes at the shortest dis- 
tance. Two (four) configurations (plus translations) were 
included in each case. The results for 4 holes on a 8 x 8 
cluster and U — 8t are, — 34.06i (four holes on a square) 
and — 34.48i (two bipolarons). Note that although the 
fourfold-polaron has also hopping terms which do not 
vanish in the infinite U limit, they are weaker than in 
the bipolaron case. These results indicate that for large 
and intermediate U no hole segregation takes place (for 
small U see below) and that the most likely configuration 
is that of separated bipolarons. 

The picture presented above is consistent with other 
analytical and numerical studies of hole-hole pairing in 
real space p5 22, 26|j2'i| at low fillings and intermediate 
to large values of U/t. As U is reduced, the size of the 
spin polarons increases and becomes elongated along the 
diagonals of the square lattice. The most likely solu- 
tions of the Hartree-Fock calculations are domain walls 
which separate antiferromagnetic regions [p8] pO|JTT| . The 
breakdown of translational symmetry associated with 
these solutions is probably real and not just an artifact 
of the Hartree Fock solution, as in the case discussed pre- 
viously. Hence, we expect a sharp transition between a 
regime of small spin polarons with an effective attrac- 
tion and striped phases at low U. Note, however, that 
the scheme presented here, based on mean field solutions 
plus corrections, is equally valid in both cases. 

Summarising, we have analyzed the dynamics of spin 
polarons and their interactions by systematically improv- 
ing the mean field approximation to the Hubbard model. 
Our scheme gives an intutive framework in which the ap- 
pearance of attraction between holes can be understood. 
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